Your thesaurus codes are: 
09.07.1, 10.08.1, 11.09.3, 11.09.4 



ASTROPHYSICS 
30.9.1994 



NGC 4526 gas, High Velocity Clouds, and Galactic Halo 
gas: the Interstellar Medium towards SN 1994D* 



D 

1 

2 
3 
4 
5 
6 
7 



L. Kingi, G. Vladilo^, K. Lipman^, K. S. de Boer'', M. Centurion^ P. Moritz^ N. A. Walton^ 

Royal Greenwich Observatory, Madingley Road, Cambridge, CB3 OEZ, United Kingdom 
Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, 34131 Trieste, Italy 
Institute of Astronomy, Madingley Road, Cambridge, CB3 OHA, United Kingdom 
Sternwarte, University of Bonn, Auf dem Hiigel 71, D-53121, Bonn, Germany 
Instituto de Astrofi'sica de Canarias, Via Lactea s/n, 38200 La Laguna, Tenerife, Spain 
Radioastronomisches Institut, University of Bonn, Auf dem Hiigel 71, D-53121, Bonn, Germany 
Royal Greenwich Observatory, Apartado 321, Santa Cruz de La Palma, 38780 Tenerife, Spain 



Received , accepted 



■ Abstract. We present spectroscopic observations of su- 
q!^ pernova 1994D in NGC 4526, an SO3 galaxy in the Virgo 
cluster 15 Mpc distant. The datasets consist of the inter- 
1— H stellar Call and Nal lines towards the supernova at high 
■4-^ spectral resolution (FWHM 6 km s~^), Ha and [Nil] ob- 
servations at lower resolution (FWHM 33 km s~^) of the 
w' nucleus of NGC 4526 and the supernova, obtained with 
the William Herschel Telescope at La Palma, and 21 cm 
04 spectra obtained with the 100 m Effelsberg Radiotelescope 

in the field of NGC 4526. 
10 The velocity of the gas in NGC 4526 determined from 
our Ha spectra is +625 km s""' at the centre (systemic 
velocity) and +880 km s~^ at the supernova position. Our 
I value of the systemic velocity is higher than the value of 
+450 km s~^ frequently quoted in the literature. 

In our high resolution spectra we detect Call and Nal 
absorption at +714 km s~^ which is produced in inter- 
C^stellar gas in NGC 4526. To our knowledge this is the first 
^ detection of interstellar absorption originating in a galaxy 
?H of early morphological type. The ratio A(Na'^)/A(Ca"'") ~ 
C/j 4 suggests an origin in cold gas at rest velocity relative to 
^ its galactic environment. The lack of multiple components 
indicates a relatively simple structure of the interstellar 
medium in the inner region of NGC 4526, at least in the 
particular line of sight to the supernova at the border of 
the nuclear ring of dust. 

We detect multi-component Call and Nal absorption 
lines in the range from +204 to +254 km s~^ which orig- 
inate in a complex of High Velocity Clouds (HVCs) lo- 



cated at a distance ^ 1 Mpc, in the surroundings of the 
Milky Way. This rare detection of HVCs in absorption 
enables the study of the properties of the gas using the 
Ca"'" and Na*^ column densities, combined with the H*^ 
column density taken from the literature at +215 km s~^ 
in the same line of sight. We find A(Na°)/A(Ca+) = 0.1 - 
0.4, in our Galaxy the signature of high velocity gas. The 
Ca"'"/H'^ and Na'^/H'^ column density ratios are extremely 
high compared to Milky Way interstellar values; the gas 
appears to have near solar abundances, very low dust con- 
tent, and a diluted ultraviolet radiation field. This is en- 
tirely consistent with Galactic fountain models, in which 
hot gas is expelled into the outer halo, and subsequently 
cools. 

At —29 km s~^, we find weak Call absorption and 
weak Hi emission. This component has properties simi- 
lar to those of the warm gas around the Sun and may 
originate in gas infalling onto the Galactic disk, perhaps 
associated with the extended complexes of Galactic halo 
gas at intermediate negative velocities which are present 
in the northern Galactic hemisphere. 

Finally, close to rest velocity, we find both warm and 
cold gas located beyond 65 pc, probably associated with 
high latitude gas at the border of Loop I. 

The total reddening of the supernova, estimated using 
the standard Milky Way gas-to-dust ratio, is E(B—Y) ~ 
0.05. 

Key words: Interstellar medium: general; Galaxy: halo; 
Galaxies: intergalactic medium, interstellar matter 
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* Based on observations taken with the William Herschel 
Telescope operated on the island of La Palma by the Royal 
Greenwich Observatory at the Spanish Observatorio del Roque 
de Los Muchachos of the Instituto de Astrofi'sica de Canarias. 



Supernovae have been studied in depth due to their impact 
on many astrophysical problems, particularly enrichment 
of the interstellar medium, neutrino physics and heavy 
element synthesis. The advent of high resolution spectro- 
graphs and efficient detectors has, over the last decade, 
allowed us to study the interstellar and intergalactic gas 
in the direction of bright supernovae in external galax- 
ies through absorption line spectroscopy with sensitivities 
previously unattainable. SN 1987A enabled the study of 
the interstellar medium in the direction of and within the 
Large Magellanic Cloud to a distance of 0.05Mpc (Vidal- 
Madjar et al. 1987, de Boer et al. 1987). SN 1993J ex- 
tended the search for interstellar and intergalactic ab- 
sorptions out to the distance of the M81 group at 3.6 
Mpc, where evidence of Galactic gas, intermediate veloc- 
ity clouds, gas within the host galaxy and tidally stripped 
material in the intergalactic space within the M81 group 
was found (Vladilo et al. 1993,1994, de Boer et al. 1993, 
Bowen et al. 1994). SN 1994D now allows us to probe gas 
out to much larger distances. SN 1994D (/ = 290.2°, h = 
70.1°) in NGC4526, an SO3 lenticular galaxy (Sandage 
1961) and a member of the Virgo cluster of galaxies at a 
distance of 15 Mpc (Pierce et al. 1994) was discovered on 
March 7, 1994 (Treffers et al. 1994) using the automated 
Leuscher Observatory Supernova Search, about one week 
before maximum brightness. Until comparatively recently, 
early type galaxies were considered to be almost free of 
dust and gas. SO3 galaxies by definition contain dust; some 
are now known to contain neutral hydrogen and molecular 
gas and have been detected at X-ray wavelengths (Roberts 
et al. 1991). The line of sight to SN 1994D is only 5° away 
from that towards SN 1991T in NGC4527 which has been 
studied using interstellar Call absorption lines (Meyer & 
Roth 1991). We therefore have the opportunity to study 
gas along a similar sight line. 

We present an analysis of high resolution observa- 
tions of interstellar Call and Nal absorption lines towards 
SN 1994D. In addition, we present lower resolution Ha 
and [Nil] observations of NGC4526 and 21 cm emission 
line profiles obtained in the field of NGC 4526. In the next 
section we describe the observations and line fitting anal- 
ysis. We discuss the results in section 3. 

2. Observations and data reduction 

2.1. Optical spectra 

The optical data were obtained using the William Herschel 
Telescope (WHT) at the Observatorio del Roque de Los 
Muchachos, La Palma, Spain. 

The high resolution data were obtained with the 
Utrecht Echelle Spectrograph (Walker & Diego 1985) at 
the Nasmyth focus of the WHT. The detector was a cooled 
Tektronix CCD with 1024 x 1024, 24 micron square pix- 
els. A spectrograph slit width of one arcsecond was used. 



D2, Di. Details of the observations are given in Table 1. 
At the time of the observations the visual magnitude was 
12.5 (Lewis 1994). We also obtained Nal spectra of two 
bright foreground stars in the field of NGC 4526 to study 
the distribution of local gas close to the line of sight to 
SN 1994D. These stars are bright (niy ~ 5), with high ro- 
tational velocity; the spectra of these stars were used to 
correct for telluric absorption in the supernova Nal spec- 
tra. 

In addition to the high resolution UES observations 
of SN 1994D, lower resolution spectra covering the Ha 
and [N ll] lines were obtained using the ISIS spectrograph 
(Carter et al. 1993), a 1200 lines per mm grating and an 
EEV CCD with 1242 x 1152, 22.5 micron square pixels, at 
the Cassegrain focus of the WHT. The one arcsecond wide 
slit, yielding a resolution of 33 km s~^ at Ha, was placed 
across the nucleus of NGC 4526 and SN 1994D. Details of 
the observations are again given in Table 1. 

Calibration frames - bias frames, fiatfields at each in- 
strumental setup and frequent arc exposures (thorium- 
argon hollow cathode lamp for the high resolution data, 
copper-neon hollow cathode lamp for the lower resolution 
data) to wavelength calibrate the data were also taken. 



Table 1. Journal of the optical observations 



Date 


Instrument 


Target 


Central 


Exposure 








A (A) 


(sees) 


09/3/94 


ISIS 


SN1994D 


6563 


1750 


15/3/94 


UES 


11 


5900 


10,800 


11 




HD 110411 


11 


180 


11 




HD 110423 


11 


360 


11 




SN1994D 


3955 


1800 


11 




11 


3953 


1800 


16/3/94 




11 


11 


7200 






11 


3961 


1800 



All optical data sets were reduced using the Starlink 
FIGARO package (Shortridge 1990). The usual steps of 
bias subtraction, fiatfielding, cosmic ray removal, wave- 
length calibration and sky subtraction were performed. 

The hollow cathode lamp frames were reduced in an 
identical manner to the object frames. The spectral reso- 
lution determined from the profile of the arc lines was for 
both the UES and the ISIS data in good agreement with 
the nominal resolving power. The data were binned to two 
bins per resolution element. 

The UES data were corrected for atmospheric absorp- 
tion in the Nal region using the template spectrum from 
the foreground stars. The individual spectra were then 
corrected to a common LSR velocity scale; each spectral 
region was then co-added and normalised. 
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Fig. 1. Interstellar lines towards SN 1994D. Components from each doublet are indicated by tickmarks. Note that there is an 
overlap in wavelength between the Nal D2, Di components. The observed data are plotted as a fine line, the synthetic profiles 
as a heavy line. 



In Table 2 we present the signal-to-noise ratio in the 
spectral regions adjacent to the absorption features. We 
also give the minimum detectable equivalent width, calcu- 
lated as Wmin = 3AAj„s(/-^. These limits are converted 
to column density upper limits using the linear part of the 
curve of growth (Spitzer 1978). 



For the ISIS data, the seeing, determined from the pro- 
file of the guide star at the time of the observations, was 
1.1 arcseconds; one dimensional spectra were created from 
the two dimensional sky subtracted spectrum by extract- 
ing in 4 spatial pixel increments (1 pixel = 0.33 arcsec- 
onds), centred on the nucleus of NGC4526 and extending 
out to the supernova position. 



Table 2. Data Characteristics 



Transitions 


S/N EW limit (3(7) 
(mA) 


N{X) limit (3(7) 
(IQi" cm-^) 


Call 


54 - 102 4.5 - 2.5 


5-3 


Nal 


75 5 


2.5 



Note:- The Call H, K data were obtained with three different 
central wavelength settings leading to varying signal-to-noise 
ratios for these lines. 



2.2. Radio observations 

The Effelsberg 100m Radiotelescope was used on April 
8, 1994 to obtain 21 cm Hi emission observations in the 
supernova field, one set centred on NGC4526 and the sec- 
ond centred ten arcminutes due north. Spatial resolution 
(half power beam width ) at 21 cm is 9.1 arcminutes and 
the main beam efficiency rjmb = 0.72. The observations 
were performed using a cooled, dual channel FET ampli- 
fier (system temperature Tgys ~ 30 K) and 1024 channel 
autocorrelator, split into two receivers with 512 channels 



of 1.29 km s~^, and a velocity range of 660 km s~^. Data 
were obtained at both positions centred on 450 km s~^, 
with integration times of 1440 seconds, and at the first 
position centred on —50 km s~^, with an integration time 
of 480 seconds. The frequency switching technique was 
used with reference frequency offsets by —3.5 MHz and 
+3.5 MHz for the velocity centres 450 km s~^ and —50 
km s""'^ respectively. All measurements were calibrated us- 
ing the brightness temperature scale calibration point S7 
(Kalberla et al. 1982) and corrected for stray radiation 
(Kalberla et al. 1980). 



3. Analysis of the spectra 

3.1. UES spectra 

The Call and Nal interstellar spectra are shown in Fig. 
1. Clearly, lines are blended; therefore profile fitting was 
performed by means of multi-component synthetic spec- 
tra. The analysis was performed using XVOIGT (Mar and 
Bailey 1994), where synthetic absorption profiles are com- 
pared with the observed profiles. Each component is repro- 
duced with a Voigt profile, using the formalism described 
originally by Stromgren (1948) - in practice, three pa- 
rameters describe the profile: radial velocity, broadening 
parameter, and column density. Combining these parame- 
ters with rest wavelength and oscillator strength (Morton 
1991) enables the formation of the intrinsic absorption 
profile which, when convolved with the instrumental pro- 
file (determined from the FWHM of unblended lines in 
the thorium-argon spectrum), can be adjusted by itera- 
tion to obtain a best fit to the observed profile. A fur- 
ther constraint is achieved by fitting both components of 
each doublet simultaneously to derive self-consistent re- 
sults. Results of the absorption line analysis are shown in 
Table 3. 

A total of seven absorption components in Call and 
five in Nal are seen, covering the radial velocity range 
—29 km s~"'^< v\sT< -|-715 km s"-'^ with good agreement 
between the two species. These break down into three dis- 
tinct groups: those with velocities in the ranges —29 km 
s""^< wisr< -5 km s"\ -1-204 km s~^< v\sr< +254 km s"\ 
and v\sT— +714 km s""'^. 



3.2. ISIS spectra 

The H a and [N ll] spectra, from the nucleus out to the 
supernova position, are shown in Fig. 2. Clearly, there 
is an increase in v\st from the nucleus outwards, with a 
corresponding change in intensity. Fits to the emission line 
profiles were performed using the ELF routines contained 
in the Starlink DIPSO package (Howarth et al. 1994). The 
fit is performed by specifying an estimate of the central 
wavelength and width of each emission line; gaussians are 
then fitted to the observed profile and then optimised by 



spectrum in an iterative manner. 

The systemic velocity of NGC 4526 as given by the 
lines formed in the nucleus extraction is +625 km s~^. 
This increases to +880 km s-i at the supernova position. 
Note that close to the nucleus of NGC 4526 the lines are 
asymmetric. 
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Fig. 2. Ha and [Nil] lines in NGC 4526. The spectra start 
at the nucleus of NGC 4526 at the bottom, incrementing in 
1.3 arcsecond steps out to the supernova position in the top 
spectrum. The rest wavelength of H a is 6562.8 A and of [Nil] 
are 6548.1 and 6583.4 Angstroms. 



3.3. Effelsberg spectra 

The 21 cm spectra of NGC 4526 are shown in Fig. 3. Gaus- 
sians were fitted to the emission line profiles using the 
CLASS package (Forveille et al. 1992). To be fitted ade- 
quately, the data required three components, the results 
of which are presented in Table 4. These components are 
in good agreement with the low velocity components seen 



Table 3. Absorption Components towards SN 1994D 







Call H & K 














Na 


I D2,Di 






Vlsi 







lo 


s, N 




Vlsi 




b 


1 „ „ AT 

log IN 


i V IN a j 


(km 




(km s"^) 




(cm- 


) 


(km 




(km s"^) 


(cm-^) 


7V(Ca+) 


-29.0 


+ 1.0 
-1.0 




11 


1 

lo 


+0 
-0 


05 

11 










< 10.40 ^ 


< 0.19 


-7.0 


+ 1.0 
-1.0 




11 


DZ 


+0 
-0 


03 
07 


-5.5 


+ 1.0 
-1.0 


D.U 


+1.0 
-1.0 




0.48 


204.5 


+ 1.0 
-1.0 


4 +^-'^ 
-4.0* 


11 


OO 


+0 
-0 


08 
05 










< 10.40 ^ 


< 0.15 


216.0 


+ 1.0 
-1.0 


5 5 +^-^ 
-3.0 


11 


OO 


+0 
-0 


03 
05 


214.0 


+2.0 
-2.0 


y .U 


+3.0 
-5.0 




0.42 


231.8 


+0.3 
-0.3 


3 5 

-0.5 


12 


1 n 

1 u 


+0 
-0 


03 
05 


234.0 


+ 1.0 
-1.0 


o . u 


+ 1.0 
-1.0 




0.12 


253.5 


+0.5 
-0.5 


2 5 +^-^ 
-1.5 


11 


55 


+0 
-0 


05 
03 


254.0 


+ 1.0 
-1.0 


2.5 


+ 1.0 
-2.0 




0.14 


714.5 


+ 1.5 
-1.5 


3 5 


11 


03 


+0 
-0 


14 
23 


713.0 


+0.5 
-0.5 


4.5 


+0.5 
-0.5 




3.55 



* Component is unresolved. 
^ 3 (7 detection limit. 



in our high resolution UES spectra. The component seen 
at —7.0 km s~^ in our UES data is split into two com- 
ponents at —8.9 km s~^ and —2.5 km s~^ in our higher 
resolution radio data. This result has been confirmed by 
observations obtained for the Leiden/Dwingeloo Hi sur- 
vey (Burton 1994). 

3.5 I < < < < < < 1 




Fig. 3. 21 cm profile towards SN 1994D 



4. Results 



We make use of kinematic arguments and the 
N(fisS')/N(Ga~^) ratio to determine where in the line of 
sight to SN 1994D the absorption components are formed. 
Specifically, calcium is more easily incorporated onto dust 
grains than sodium; in galactic disks, where gas is gener- 
ally quiescent, #(Na°)/#(Ca+) > 1, but if the dust grains 
are disrupted, either by shock waves or cloud collisions, 
Ca"'" is released into the gas phase and N(fisS')/N(Ga~^) 
< 1. This latter condition is frequently found in clouds 
moving at high velocity ( > 20 km s~^) relative to the 
local rest velocity, both within a galactic disk and distant 
from it (Berlin et al. 1993, and references therein). It is 
also possible to find a low N(fia°)/N(Ga~^) ratio in gas 
at low velocity if this gas is warm as Na*^ is preferentially 
removed by collisional ionisation. 



Table 4. Hi gas towards SN1994d 



Vlsi 


FWHM 


N{R") 


T 

-1- max 


(km s-i) 


(km s-i) 


(10i« cm-2) 


(K) 


-8.9 ±0.1 


18.7 ±0.2 


115.5 ±1.7 


7540 


-2.5 ±0.1 


4.7 ±0.4 


7.3 ±0.7 


470 


-29.2 ±0.9 


17.3±1.7 


14.3 ±1.6 


6520 



Tmax- upper limit assuming pure thermal broadening of the 
HI emission 



Atomic hydrogen originating in NGC 4526 was generally 
undetected at 21 cm (Kumar & Thonnard 1983, Giova- 
nardi et al. 1983), with the exception of the observations 
obtained by Davies & Lewis (1973) which yield a central 
velocity of +450 km s~^, a value frequently quoted in the 
literature. We do not detect atomic hydrogen near this 
velocity in our Effelsberg data, or in our low resolution 
H a or [N ll] optical spectra. Our H a spectra yield radial 
velocities ranging from +625 km s~^ at the nucleus of 
NGC 4526, up to +880 km s"-'^ at the supernova position. 
A similar value for the central velocity, t;hei = +602 km 
s~^ was found by Huchra et al. (1983) and Huchra (1987), 
from optical spectroscopy as part of the CfA redshift sur- 
vey (t^hei = fisr+ 3.8 km s""'^). A higher range of radial 
velocities for the interstellar gas in NGC 4526 is also sup- 
ported by the CO emission measured by Sage & Wrobel 
(1989) which peaks at t;hei = +743 km s"-'^ and extends 
from +320 km s'^^ vx^^ <, +1060 km s'^. 

The kinematic evidence suggests that the component 
present in the high resolution spectra with visr = 714 
km s""'^ is produced in NGC 4526; the components at visr 
300 km s~^ are however unlikely to originate in that 
galaxy. 

The ratio #(Na°)/#(Ca+) ~ 4 (Table 3) for this com- 
ponent is a typical value for interstellar gas close to rest 
velocity in our Galaxy (Siluk & Silk 1974). This, coupled 
with the good radial velocity agreement suggests that the 
visr = +714 km s""'^ component originates in rest velocity 
gas in NGC 4526. 

NGC 4526 is a lenticular galaxy classified as SO3 owing 
to the presence of a nuclear ring of dust (Sandage 1961). 
Comparison of blue and near-infrared CCD images shows 
that the dust ring has an angular extent of ~ 31" x 6" 
(Veron-Cetty & Veron 1988). The detection of molecular 
gas (Sage & Wrobel 1989) and the fact that NGC 4526 
is a bright IRAS galaxy (Soifer et al. 1989) are probably 
due to the presence of this dust. The line of sight towards 
SN 1994D, at 7" North and 9" West of the nucleus (Tref- 
fers et al. 1994), lies a few arcseconds from the NW border 
of the ring; the NGC 4526 interstellar gas at +714 km s""'^ 
is probably associated with the outskirts of the dust ring. 

By adopting the average #(Na'^)/#(H) ratio of the 
Milky Way interstellar gas (Ferlet et al. 1985) we estimate 
#(H) ~ 7.6 X 10^^ cm~^ for the total gas column density 
in the line of sight within NGC 4526. This value is rela- 
tively low for a line of sight which intersects a galaxy close 
to its nucleus, and is consistent with the low gas content 
expected for SO galaxies, which are often (~ 80% of the 
cases) undetected in Hi emission (Roberts et al. 1991). 
Moreover, the absence of multiple components suggests 
that the structure of the interstellar medium of the host 
galaxy is relatively simple, at least in the line of sight to- 
wards SN 1994D. Clearly care must be taken in extending 
these results to the galaxy as a whole and, more generally. 



nents and total column density originating in NGC 4526 
may have been much higher if the line of sight had inter- 
sected the nuclear ring along its maximum depth rather 
than at the edge. 

4.2. High Velocity Clouds 

The multiple absorption complex detected in Call and 
Nal at +204 <, vi^r ^ +254 km s"! has #(Na°)/#(Ca+) 
between 0.1 and 0.4, typical of high velocity gas (Siluk 
& Silk 1974). The 21 cm counterpart of this gas is not 
detected in the Effelsberg data at a 3(7 detection limit 
#(H°)um ~ 7 X lO^s cm-2. Kumar & Thonnard (1983) 
however, measured weak emission around ~ +215 km s~^ 
with N{E°) ~ 2 X 10^® cm'^ by pointing the Arecibo 305 
m radiotelescope (FHWM = 3.3') in a strip of ~ 10' x 3.3' 
along NGC 4526. The constant radial velocity of this fea- 
ture along NGC 4526 indicates that the gas is not asso- 
ciated with that galaxy. As discussed previously, the ve- 
locities of the multiple absorption complex being ^ 300 
km s~^ indicate it is unlikely they are physically associ- 
ated with NGC 4526. Furthermore, since this cloud com- 
plex has a mean velocity of ~ —400 km s""'^ with respect 
to the systemic velocity of NGC 4526 and of ~ +160 km 
s""'^ with respect to that of the Milky Way (vqc = ^isr 
+ 220 sin / cos h km s""'^), it is more probably related to 
the surroundings of our Galaxy than to the environment 
of NGC 4526. The angular dimension of the order of 10' 
argues against a location in intergalactic space at a dis- 
tance beyond 1 Mpc, since this would imply that the gas 
has linear scales of 3 kpc, with a column density of only 
2 X 10^^ cm~^, i.e. the gas would have a local density of 

2.2 xlO"'* atoms cm~^ - even lower than that of the 
hot phase of the interstellar space where species like Na*^ 
cannot survive. A distance in the range from 1 to 10 kpc, 
corresponding to a linear size of 3-30 pc, would give a 
local density of 0.22-0.02 atoms cm~^, within the range 
of normal interstellar cloud sizes and densities. 

The cloud complex at +204 ^ visr ^ +254 km s~^ 
cannot be gas corotating in the Milky Way halo, since in 
this case it would have |t;isr| ^ 5 km s""'^ owing to the 
relatively small angle between the direction of the super- 
nova (/=290°) and that of the solar motion around the 
Galactic centre (see Fig. 4). The complex is not at rest 
with respect to our Galaxy — if it were at rest it would 
have Visr = +70 km s""'^ as a result of the projection of 
the solar motion along the line of sight. The gas complex 
is most likely related to the phenomenon of High Veloc- 
ity Clouds (HVCs). Although the origin and location of 
HVCs is not yet fully understood, it is clear that many 
lie in the surroundings or even in the halo of our Galaxy 
(Wakker & van Woerden 1991). Only a small number of 
HVCs were detected by Hulsbosch k Wakker (1988) in a 
20°x20° field centered on NGC 4526 at a detection limit 
of #(H'^) ^ 2 X 10^^ cm~^; all were at positive velocities. 



general span a very wide range of negative and positive 
velocities, the relatively small velocity difference between 
the gas at +204 ^ visr ^ +254 km s"^ and the HVCs in 
the field suggests that we are dealing with HVCs unde- 
tected in the survey of Hulsbosch & Wakker (1988). Our 
identification of these components as HVCs is supported 
by the detection of HVCs at remarkably similar veloci- 
ties (215 - 263 km s'^) towards SN 1991T in NGC4527 
(Meyer & Roth 1991), 5° from the SN 1994D sightline - a 
distance compatible with the size of known HVCs. Alter- 
natively, it is possible that the line of sight to SN 1994D 
intersects a small Very High Velocity Cloud; several were 
found by Hulsbosch & Wakker (1988) in the northern sky 
with velocities in excess of 200 km s""'^. 

Our Call and Nal detections, together with the Hi 
detection of Kumar & Thonnard (1983) provide the rare 
opportunity to compare absorption and emission col- 
umn densities of HVCs. The ratios N(Ca+)/N(R°) and 
N(fisS')/N(li°) of the summed components in the range 
from -|-204 to -|-254 km s""'^ are extremely high with 
respect to those found in the ISM of our Galaxy. We 
find log N(Ca+)/N(R°) = -6.0 dex, while the highest 
value from the compilation of interstellar column densi- 
ties of Stokes (1978) is —8.1 dex for disk gas; our value 
is also higher than that found in more recent work (Al- 
bert et al. 1993), concentrating on interstellar absorp- 
tion towards halo stars. Even considering the lowest Ca"'" 
ionisation fraction found by Phillips et al. (1984), i.e. 
#(Ca++)/#(Ca+) ~ 1, this implies log N(Ca)/N(R°) 
= —5.7 dex for the fraction of the total calcium, a value 
remarkably close to the solar abundance, [Ca/H] = —5.64 
dex (Anders & Grevesse 1989). 

For sodium, we find #(Na°)/#(H°) = -6.8 dex. In 
comparison, the average Galactic value is —8.3 dex (Stokes 
1978) and the solar abundance is —5.71 dex. 

To bring together these results, we explore three sce- 
narios. Firstly, if we accept the calcium abundance per se, 
then either the gas is of solar metallicity with minimal 
depletion, or, if there is depletion, the gas must be even 
more enriched. This proposal is consistent with the below- 
solar column density of sodium only if the majority of it 
is present in an unobserved ionisation state such as Na"'" 
- an expected result due to the low ionisation potential 
of Na*^, = 5.14eV. In fact, it is surprising that we see so 
much Na*^ compared to the mean value in the ISM. This 
suggests we are seeing an ionisation balance for sodium 
different to that in the local ISM, which may be a symp- 
tom of a diluted UV radiation field which is expected far 
from the Galactic plane. 

We must recognise however that the detailed ionisa- 
tion balance of this gas is largely unknown, and to at- 
tribute our results to any single ionisation mechanism is 
extremely precarious. This leads us to the second possibil- 
ity: that the gas is, to a significant extent, ionised. In this 
case, the total hydrogen column density would be substan- 



the metallicity of the gas would be lower than the near- 
solar values suggested in the previous scenario. Enhanced 
ionisation however fails to explain the excess of neutral 
sodium compared to that seen in the Galactic ISM. 

Lastly, a significant fraction of the hydrogen could 
be in molecular form, which would again imply that the 
metallicity of this gas is below solar. However, such an 
explanation is deemed unlikely since there is no known 
association between molecular gas and HVCs. 

The most likely of these situations appears to be the 
first: that is, the gas has near solar abundances, with min- 
imal depletion onto dust. This is in agreement with the 
proposal of Lipman & Pettini (1994), based on observa- 
tions of Till towards halo stars, that gas ^ Ikpc from the 
Galactic plane contains little dust and has essentially zero 
depletion. 

As mentioned previously, absorption at similiar veloci- 
ties has been detected towards SN 1991T, a sightline 5° 
away (Meyer & Roth 1991). These components are an 
order of magnitude weaker than those we have observed 
towards SN 1994D; since no absorption was detected by 
the same authors towards 3C273, at 1° further iiom the 
SN 1994D sightline, we conclude that if one HVC complex 
is responsible for the Call absorption in this field, the 
SN 1991T and 3C273 sightlines intersect its low column 
density periphery. 

To date there is little information available regard- 
ing the abundances found in HVCs. Danly et al. (1993), 
studied the absorption lines of Cll, Sill and Ol and re- 
port a lower limit for the metallicity of HVC complex M 
at TiHVC > —1.0 to —1.4 dex of solar, but predict that 
the gas is considerably more enriched than this. Bowen 
& Blades (1993) detected two HVC components towards 
Mrk205 using Mgll absorption lines; they found that one 
had TiHVC = —2.0 dex of solar (though this may represent 
a lower limit due to the possible presence of unresolved 
high column density components) and for the other they 
derived a lower limit of Tihvc = +0.4 ± 0.5 dex of so- 
lar. Robertson et al. (1991), Songaila (1981) and West et 
al. (1985) detected Call absorption arising from known 
HVCs; the N(Ca+)/N(H°) ratios yield depletions of -1.2 
to —2.7 dex. These values assume an insignificant amount 
of calcium is present as Ca"'""'" and therefore overestimate 
the depletion. If we assume N{Ga++)l N{Ga+) ~ 4, the 
mean value in our Galaxy, the depletions are —0.5 to —2.0 
dex. 

With such limited information it is impossible to accu- 
rately account for the presence of elements in unobserved 
ionisation states; nevertheless it appears that an apprecia- 
ble fraction of HVCs may have near-solar abundances, and 
suffer from little depletion. Observations of different ele- 
ments and ionisation species will be necessary in order to 
gain greater insight into the ionisation balance and dust 
depletion in these clouds more fully, and hence confirm 
this hypothesis. 



with Galactic fountain models (Bregman 1980, Spitzer 
1990) in which hot gas is expelled high into the halo and 
subsequently cools. Such an energetic process would re- 
turn the refractory elements from the dust grains back to 
the gas phase, explaining the low depletion of calcium, 
and could also explain the high overall metallicity and 
anomalous velocities of the HVC gas. The typical turning 
height of the fountain models is a few kiloparsecs; as we 
only see clouds moving away from the Galactic plane and 
not inf ailing, this is consistent with our estimation of the 
distance to the clouds being in the range 1 to 10 kpc. 

4.3. Milky Way interstellar gas 

Fig. 4 shows the expected range of function of dis- 

tance from the Sun in the direction of SN 1994D assuming 
circular rotation and the Galactic rotation curve of Brand 
and Blitz (1991). Clearly, nearly all the components seen 
in our high resolution spectra have velocities which can- 
not be explained by the Galactic rotation curve. The two 
components at v\st = —7 km s""'^ and —29 km s""'^ are most 
likely produced in Milky Way or Halo gas. These compo- 
nents are not visible in the spectrum of the foreground 
stars HD 110411 or HD 110423, setting a lower limit on 
the distance of 65 pc to absorbing clouds. 
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Fig. 4. Vlsr as a function of distance from the Sun in the 
direction of SN 1994D assuming the rotation curve of Brand 
& Blitz (1991) and circular rotation. The effect of rotation 
parameters deviating from the reference curve can be found in 
Kaelble et al. (1985) 



The component at —29 km s~^ is detected in Call but 
not in Nal. It is also detected as an asymmetry on the 
blue wing of the main emission at 21 cm. The width of 
this 21 cm component indicates that the gas is warm or 
highly turbulent (Table 4). The low values of the neutral 
hydrogen and neutral sodium column densities - #(H'^) 
~ 1.4x10^^ cm-2 and #(Na°) < 2.5 x 10^° cm'^ - suggest 
that we are dealing with warm gas if the correlations be- 
tween physical parameters and column densities found by 
Vladilo & Centurion (1994) hold for this component. The 
ratio #(Na°)/#(Ca+) < 0.19 is very similar to that found 
by Berlin et al. (1993) for the warm cloud in the vicinity of 
the Sun. The ratio log #(Ca+)/#(H°) = -8.0 dex is high 
compared to the typical disk values which are in the range 
from —9.7 to —8.1 (see Stokes 1978), but well within the 
range from —8.6 to —6.5 found in the Galactic Halo (Al- 
bert et al. 1993). Since the interstellar Ca+^/C a"'" number 
ratio ranges from 1 to 10 with a mean value of 4 (Phillips 
et al. 1984), we find log #(Ca)/#(H°) = -7.4 (t°;^) dex. 
This implies that calcium is depleted by —1.8 (Iq 5) '^^^ 
with respect to the solar abundance value log [Ca/H] = 
—5.64 dex (Anders & Grevesse 1989). This value for the 
depletion is relatively low for calcium and is close to the 
depletion found by Berlin et al. (1993) for the local inter- 
stellar cloud. The negative radial velocity of the —29 km 
s~^ component cannot be explained in terms of Galactic 
differential rotation either assuming co-rotation or slower 
rotation at large distances from the Galactic plane (see 
Kaelble et al. 1985). The gas at —29 km s~^ has therefore 
an intrinsic velocity component directed towards the disk 
of our Galaxy; it is probably related to the phenomenon of 
intermediate negative velocities found in the general direc- 
tion of the North Galactic Pole (Wesselius & Fejes 1973). 
This gas is likely to be part of the infalling structure lo- 
cated between 250 and 1700 pc above the disk proposed 
by Danly (1989). 

4.3.2. Local gas 

The low velocity of the absorption component at v\st 
= —7 km s""'^ suggests that it is produced in nearby 
gas. The ratio #(Na'^)/#(Ca"'") = 0.5 however is low 
for gas close to rest LSR velocity, which has typically 
#(Na°)/#(Ca+) > 1 (Siluk k Silk 1974). The low value of 
the #(Na'^)/#(Ca"'") ratio in this low velocity component 
can be explained if the gas is warm, as is the case for the 
local interstellar medium in the solar vicinity (Berlin el 
al. 1993). The 21 cm emission profiles indicate that the —7 
km s~^ absorption is in fact a blend of two components, 
one narrow component originating in cold gas and one 
broader component which may originate in warm gas (see 
Table 4). The #(Ca+)/#(H°) and #(Na°)/#(H°) ratios 
(sum of the two components) are both in the range typical 
of Galactic clouds, as expected for local gas. The —7 km 



it is not seen in front of the foreground star HD 110411 
(/=295°,&=+73°) at 65 pc distant, the gas at -7 km s'^ 
lies at least 60 pc from the Galactic plane, probably be- 
yond the edge of the Local Bubble (Cox & Reynolds 1987). 
The gas may be associated with the northern border of 
Loop I (Berkhuijsen et al. 1971) which is intersected by the 
supernova line of sight. The tangential distance of Loop 
I is 70 ± 40 pc (Berkhuijsen 1973) consistent with the 
lower limit for the distance suggested by our data. The 
line of sight to SN 1994D also intersects the outskirts of 
the cloud IREC 421 (Desert et al. 1988), which is centered 
at /=290.8°, &=72.7°and has an extension of 3.3 square de- 
grees. Most of the objects in the catalog of Desert et al. 
(1988) are probably high latitude molecular clouds and 
the component at —2.5 km s~^ could represent cold gas 
associated with IREC 421. 

4.4- Reddening towards SN 1994D 

The £'(B— V) reddening of the local gas and of the 
HVCs can be estimated from the neutral hydrogen col- 
umn densities by using the the average Galactic ratio 
N{R°)/E{B-Y) = 4.8 X 10^1 cm'^ (± 50%) of Bohlin 
et al. (1978). We derive E(B-Y) = 0.029 and 4.2 xlO"'' 
for the local gas and for the HVCs respectively. As far 
as the gas in NGC 4526 is concerned we are forced to 
use the indirect determination of the total hydrogen col- 
umn density (Sect. 3.1) in conjunction with the relation 
N(Rtot)/E(B-Y) = 5.8 x lO^i cm'^ (Bohlin et al. 1978). 
This gives £'(B— V) = 0.013 for the reddening due to the 
-|-714 km s~^component. The total reddening in front of 
the supernova is therefore £'(B— V) = 0.046, with the ma- 
jor source of uncertainty being the NGC 4526 interstellar 
contribution. 

5. Summary 

The line of sight to SN 1994D in NGC 4526 samples several 
interstellar media which can be summarised as follows:- 

- the interstellar gas in NGC 4526. #(Na°)/#(Ca+) is 
typical of gas at rest in the disk of a galaxy. It has 
a simple velocity structure, and is most likely associ- 
ated with the dust ring which surrounds the nucleus of 
NGC 4526. The estimated total column density is low, 
as is expected in SO galaxies. This is the first positive 
detection of interstellar absorption in a galaxy of early 
morphological type. 

— High Velocity Clouds in the range 1-10 kpc from 
us, with a linear size of ~ 3-30 pc. The depletion of 
calcium and sodium in these clouds is less than that 
normally found in the ISM of our Galaxy, they may be 
experiencing a diluted radiation field, and they appear 
to have very low dust content. This is consistent with 
Galactic fountain models, where refractory elements 
are returned to the gas phase. 



negative velocities, possibly associated with a struc- 
ture located between 250 pc and 1700 pc above the 
Galactic disk. 

— Two low velocity components formed in nearby gas. 
The N(fisS')/N(Ga~^) ratio of the combined compo- 
nents indicates that the bulk of the gas is warm rather 
than shocked gas; it is located at least 65 pc from 
us. One component is formed in cold gas, and may 
be formed in a high latitude molecular cloud. 
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